Mon. Not. R. Astron. Soc. 000, 000-000 (0000) Printed 3 March 2013 (MN MfcjX style file v2.2) 



Tracers of Discs and Winds around Intermediate and High 
Mass Young Stellar Objects 



(N 

o : 



S.L. Lumsden 1 , H.E. Wheelwright 1 ' 2 , M.G. Hoare 1 R.D. Oudmaijer 1 and J.E. Drew 

1 Department of Physics and Astronomy, University of Leeds, Leeds LS2 9JT, UK 

2 Max Planck Institute for Radio Astronomy, 53121 Bonn, Germany 

3 Department of Physical Sciences, University of Hertfordshire, Hatfield, Hertfordshire, ALIO 9AB, UK 
Email - s.l.lumsden@leeds.ac.uk 



3 March 2013 



43 

9 L,: 
6 



> 

(N 
(N 

iri 
O 

(N 



ABSTRACT 

We present a study of the kinematical properties of a small sample of nearby near- 
infrared bright massive and intermediate mass young stellar objects using emis- 
sion lines sensitive to discs and winds. We show for the first time that the broad 
(~ 500kms _1 ) symmetric line wings on the HI Brackett series lines are due to Stark 
broadening or electron scattering, rather than pure Doppler broadening due to high 
speed motion. The results are consistent with the presence of a very dense circum- 
stellar environment. In addition, many of these lines show evidence for weak line 
self-absorption, suggestive of a wind or disc-wind origin for that part of the absorb- 
ing material. The weakness of the self-absorption suggests a large opening angle for 
such an outflow. We also study the fluorescent 1.688/im Fell line, which is sensitive 
to dense material. We fitted a Keplerian disc model to this line, and find reasonable 
fits in all bar one case, in agreement with previous finding for classical Be stars that 
fluorescent iron transitions are reasonable disc tracers. Overall the picture is one in 
which these stars still have accretion discs, with a very dense inner circumstellar envi- 
ronment which may be tracing either the inner regions of a disc, or of a stellar wind, 
and in which ionised outflow is also present. The similarity with lower mass stars is 
striking, suggesting that at least in this mass range they form in a similar fashion. 

Key words: accretion discs; line: profiles; stars: wind, outflows; stars: pre-main- 
sequence; stars: formation 



1 INTRODUCTION 

There are well recognised and long-standing difficulties in 
describing high-mass star formation as a simple scaled-up 
version of the low-mass case (eg Kahn 1974, Shu, Adams 
& Lizano 1987). Radiation pressure is sufficient to prevent 
further accretion once fusion commences in the stellar core, 
preventing the formation of a star much more massive than 
20 Mq in the simple low mass theory. This led to alternative 
models of high-mass star formation being suggested (e.g. 
Bonnell, Vine & Bate 2004, Stahler, Palla & Ho 2000, Mc- 
Kee & Tan 2003) which can generally be subdivided into 
two groups: those that overcome the mass limit by form- 
ing massive stars through competitive accretion in dense 
cluster cores and those that look for a working version of 
the scaled-up low-mass star formation paradigm, essentially 
overcoming the pressure from the star through self-shielded 
disc accretion. Dynamical models of disc accretion around 
massive stars are now being produced which appear to over- 
come the mass limit (eg Krumholz et al. 2009 who model the 
formation of a 40 Mq star, and Kuiper et al. 2010 who find 



radiation pressure is not an efficient means of preventing 
disc accretion even in very massive stars). Such models are 
therefore currently favoured as the most likely. This agrees 
with limited observational data, where infrared interferome- 
try has revealed disk- like structures (eg. Kraus et al. 2010). 
However we still have little direct evidence for accretion in 
a large sample of young massive stars. 

The obvious alternative method to apply to a large sam- 
ple is to look for kinematical evidence of a disc. The 2.3/im 
CO first overtone rotation-vibration emission is believed to 
mostly arise in an accretion disc in young stars (eg, Carr 
1989, Chandler, Carlstrom & Scoville 1995, Bik & Thi 2004, 
Blum et al. 2004, Wheelwright et al. 2010, Davies et al. 
2010). The physical conditions required for its origin are 
characteristic of the properties of a disc that is self-shielded 
(temperature in the range 1000-5000K, density greater than 
10 10 cm -3 ), since it allows for the presence of molecular gas 
relatively near the exciting star. Generally, the spectra of 
objects studied in this way have been better fitted by a 
disc model than other alternatives. However not all mas- 
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sive young stellar objects show CO bandhead emission (eg 
Porter, Drew & Lumsden 1998). There are many reasons 
why this might be, but the most plausible are that the ac- 
cretion rate may simply be too low to give detectable emis- 
sion, or, there may be too little column density to create 
the self-shielding condition required in the disc to prevent 
the intense UV flux from more massive stars reaching the 
mid-plane and destroying the CO. The CO bandhead is also 
relatively complex to interpret without high signal-to-noise 
data (eg Wheelwright et al) because of the blending of the 
multiple rotational transitions that make up the v = 2 — 
vibrational transition. 

Another complicating factor in high mass stars is that 
other emission mechanisms are also likely to be present. The 
most obvious of these is a stellar wind which can clearly be 
seen in at least some massive young stellar objects (eg Drew, 
Bunn & Hoare 1993, Bunn, Drew & Hoare 1995, the models 
presented by Sim, Drew & Long 2005, and see Section 3.2). 
A powerful stellar wind is a fairly natural consequence of a 
high accretion rate in any hot star system. Winds in young 
stars are usually traced by recombination lines of common 
elements, particularly hydrogen and helium (eg Simon et 
al. 1983 and Persson et al. 1984 at moderate spectral res- 
olution, and Bunn, Drew and Hoare 1995 at high enough 
resolution to fully sample the line profiles). These lines can 
also be excited in nebular gas and in a disc as well. An ideal 
study should combine tracers sensitive to the winds as well 
as those that sample density and temperature and ionisation 
ranges more typical of an accretion disc in order to clearly 
distinguish the origin of the emission. Therefore we need a 
transition that is only strongly excited at the higher densi- 
ties seen in circumstellar discs around young stellar objects 
(eg n > 10 6 cm -3 ). 

Fortunately, there are other similar well studied situa- 
tions from which we can draw likely candidates. For exam- 
ple, spectra of both active galactic nuclei (eg Baldwin et al 
2004) and classical Be stars (Zorec et al. 2007) show fluo- 
rescent Fe + emission consistent with an origin in a disc. It 
is not seen in ordinary nebular gas (eg as in an HII region: 
see, for example, Lumsden & Puxley 1996). Fell emission 
straddles the boundary between where CO emission might 
still arise (neutral material) and ionised material since its 
ionisation potential is only 7.9eV, with Fe 2+ having an ioni- 
sation potential of 16.2eV. Carciofi & Bjorkman (2006) have 
successfully modelled the origin of this emission in Be star 
disks from an origin in the shielded material near the disk 
mid-plane. Accretion discs around young stars are likely to 
have additional opacity due to dust, but a similar partially 
ionised zone will exist within them. 

Fluorescent excitation typically occurs when an 
atom/ion is irradiated by either UV continuum or Lyman 
series photons. The excited atom/ion cascades back down 
to the ground state, leading to transitions that would not 
typically be seen from the standard recombination process. 
The optical Fell lines have been well studied and modelled 
in both active galaxies and Be stars (eg Baldwin et al. 2004; 
Carciofi & Bjorkman 2006) but there are also infrared tran- 
sitions that are more suitable for studying heavily embedded 
young stellar objects. In particular there are two strong lines 
near 1.6^m (eg see the spectrum of 77 Carina in Hamann et 
al. 1994), of which the transition at 1.6877808/im is both in a 
region of good atmospheric transmission, and near enough 



in wavelength to the HI Br 11 line that both can be com- 
pared without the need to worry about differential extinc- 
tion. The 1.6877808/xm c 4 F 9 -z 4 F 9 transition arises from an 
upper level 6.2eV above the ground state (Johansson 1978) 
and is likely to be pumped by Lya in situations just such 
as we are seeking (the dense ionised/neutral boundary) as 
discussed by Johansson & Letokhov (2007). Previous obser- 
vations (Porter et al. 1998) indicate this particular line may 
be relatively common in massive young stellar objects, and 
more prevalent there than in lower mass objects (Hamann 
& Persson 1992). 

The aim of this paper therefore is straightforward. We 
have observed a small sample of eight intermediate and high 
mass young stellar objects with sufficient spectral resolution 
to resolve the kinematics of the line emission in both Fell 
and HI recombination lines for the first time. The bulk of 
this sample have near and /or mid-infrared interferometric 
data available, allowing plausible limits to be placed on the 
inclination of the system and the size of any infrared emit- 
ting region in a disc. The comparison of both Fell and HI 
tracers is ideal for testing whether the Fell samples a dif- 
ferent emission region from the HI recombination lines. The 
sample is therefore ideal for testing the value of the Fell line 
as a disc tracer, as well as expanding on the use of HI lines 
as studied previously (eg Bunn et al 1995) to those present 
in the H band. Our work presents a velocity resolved study 
of the higher Brackett series members 11-4 and 12-4 for the 
first time. 

In section 2 we discuss the sample and the observations 
made, in Section 3 we discuss the observational results to- 
gether with some simple model fits to our data and in Sec- 
tion 4 we sum up the evidence as to which features might be 
the best kinematical tracers of discs or winds around young 
massive stars. 



2 OBSERVATIONS 

The data were acquired on the nights of 24 and 25 June 2000 
and 7 and 8 July 2001 using the facility near infrared spec- 
trometer CGS4 on UKIRT. We used the echelle grating and 
300mm camera with a one pixel wide slit. Due to anamor- 
phic distortion within the instrument the spatial scale along 
and across the slit vary. The effective size of a pixel with 
this setup is 0.4 arcseconds across the slit, and 0.91 arcsec- 
onds along it. The measured spectral resolution from an arc 
calibration line was ^krns" 1 . The slit was oriented at the 
default north-south position angle. We observed two sepa- 
rate wavelength ranges for all of our targets, from approxi- 
mately 1. 636-1. 648/im in order to detect the 1.6411725/im 
Brl2 and 1.6440018/im a 4 F 9 -a 4 D 7 [Fell] lines, and 1.679- 
1.69lAim to detect the 1.6811164/mi Brll and 1.6877808^m 
c 4 F 9 -z 4 F 9 Fell lines. We also observed the 2.159-2.173/im 
range in order to measure 2.166127/im Br7 for five of the tar- 
gets. All wavelengths are in vacuum, and hereafter all wave- 
lengths are truncated to three decimal places. All eight of 
the main targets are known massive YSOs or Herbig Ae/Be 
stars. Some basic properties of these are given in Table [1] 
Where possible we have derived inclinations from near or 
mid-infrared interferometry. For GL 490 we instead use re- 
sults from mm CO interferometric measurements, and for 
S106 IRS we rely on models of the velocity structure. Where 
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KJ UJCCL 


Distance (kpcj 


Spectral type 


Luminosity ( L0) 


Ay (mag) 


Inclination 


V LSR (kms 1 ) 


BD+40°4124/V1685 Cyg 




B3V 2 




3.1 3 


40 - 50° 4 


~8 5 


GL 490 






2000 6 


42 7 


~ 30° 8 


_ 14 9 


M17SW IRS1 


1.3/2. 1 14 


09V 15 




17 16 


76° 15 


~17 17 


M8E 


1.25 10 




20000 11 


39 s 


~ 20° 12 


~ ll 13 


MWC 297 


0.25 18 


B1.5V 18 




gl8 


5° 19 , < 40° 20 , 15° 21 


2 - 20 18 


MWC 349A 


1.7 22 




60000 23 


9 23 


> 70° 24 


10 25 


S106 IRS 


0.5 - 1.7 26 


07.5 V 26 




13.6 27 


~ 75° 28 


~5 29 


VV Ser 


0.2 - 0.6 30 


B6V 2 /B9V 31 




~3 31 


>80° 32 


~ 10 33 



Table 1. literature. 1: van den Anckcr, de Winter & Tjin A Djie (1998); 2: Hernandez et al. (2004); 3: Oudmaijer, Busfield & Drew 
(1997), van den Ancker, Wesselius & Tielens (2000); 4: Eisner et al. (2003, 2004); 5: Fuente et al. (1990) from NH 3 ; 6: Schreyer et al. 
(2002); 7: Willner et al. (1982); 8: Schreyer et al. (2006); 9: Waterlout & Brand (1989) from molecular gas (CO); 10: Prisinzano et al. 
(2005) assuming it is at the same distance as the resolved stellar cluster found in M8; 11: Thronson, Loewenstein & Stokes (1979); 12: 
Linz et al. (2009); 13: Bunn et al. (1995) from the star and Simon et al. (1984) from CO; 14: Estimates vary from 1.3kpc (Hanson, 
Howarth & Conti 1997) using spectral parallaxes to 2.1kpc (eg Hoffmeister et al. 2008) using photometric methods - we adopt the latter; 
15: Follert et al. (2010); 16: Porter et al. (1998); 17: Bunn et al. (1995) from the star, Lada (1976) from CO; 18: Drew et al. (1997); 
19: Alonso-Albi et al. (2009); 20: Acke et al. (2008); 21: Malbet et al. (2007); 22: Meyer, Nordsieck & Hoffman (2002); 23: Cohen et al. 
(1985); 24: Danchi, Tuthill & Monier (2001); 25: Gordon et al. (2001); 26: Schneider et al. (2007), we adopt the further distance, which 
their analysis shows results in the given spectral type; 27: van den Anckcr ct al. (2000); 28: Hippclcin & Mucnch (1981), Solf & Carsenty 
(1982); 29: Bunn et al. (1995); 30: Montesinos et al. (2009); 31: Pontoppidan et al. (2007); 32: Eisner et al. (2003); 33: Corcoran & Ray 
(1997). 



required the tabulated luminosities have been modified from 
the original references to take account of revised distances. 
For sources where an accurate spectral type can be inferred 
we have given that. For those where only a bolometric lu- 
minosity exists (generally the most heavily embedded) we 
have quoted that instead. All of these sources were chosen 
to be bright enough at H to be observable with CGS4. In 
addition, we observed the planetary nebula NGC 7027 and 
the yellow hypergiant IRC +10420 for calibration purposes. 
NGC 7027 provided a check of our wavelength calibration, 
since its velocity is known to a greater degree of accuracy 
than most of our target stars. IRC +10420 was observed as 
it is a known Fell emitter in the optical, and hence pro- 
vides a check of our line identification. The sample is not in 
any sense complete or statistically representative of massive 
YSOs as a whole. 

All the sources are sufficiently compact that they could 
be observed by nodding the telescope such that the sources 
remained on the slit (the measured beam separation is 22 
arcseconds). We observed standard stars to provide crude 
flux calibration and correct for atmospheric absorption fea- 
tures in addition to observations of the targets. The stan- 
dard stars were chosen to have a range of spectral type from 
G6V through to 09V. This provided enough range to enable 
us to correct for the observed Brackett line absorption in 
the standards themselves. Cooler stars were not used since 
they have many other spectral features in the wavelength 
ranges observed, and in practice we used the hottest stan- 
dard available for any target. The Brackett line absorption in 
the standard was interpolated over to remove it. This gives 
good though not perfect correction: Brll in particular may 
be affected by poor correction in the standard star since the 
line lies near the edge of the spectrum (cf Figure [1} . There 
are some hints that the Brll and Brl2 line profiles differ in 
some objects (eg BD+40°4124, Figure [TJ, but the difficulty 
in correcting adequately for the instrumental response and 
interpolation of Brll in the standard star makes it impossi- 
ble to be sure these are real. We therefore do not use Brll 
in our analysis. For Brl2 and Br7 the interpolation process 



is far more secure, and is certainly adequate to define weak 
line wings. 

A standard reduction procedure was adopted for all of 
the data discussed here. The projected image of the slit on 
the array is curved. This was corrected for using the trace of 
a bright arc line as a model. This procedure was checked for 
reliability using fainter night sky lines present on each indi- 
vidual exposure. Each frame was also checked to ensure that 
no velocity shifts were apparent, again using the night sky 
lines. Where such shifts were seen the data were realigned 
using the centroid of the night sky line. 

Both of the beams were extracted from the data after 
correction for spectral distortion with wavelength. The typ- 
ical extraction is 6 spatial pixels (or 4.5 arcseconds) to en- 
compass the full image of the spectrum as seen on the array. 
The telescope was slightly defocused to correct for problems 
caused by the order sorting filter present in CGS4, which 
utilises a circular variable filter (CVF). This imposes a wave- 
length dependent transmission pattern on the spectra that 
contains both a high frequency component (approximately 
one resolution element) and a lower frequency component 
where the CVF response declines away from the centre of 
the observed wavelength range, giving lower transmission 
and signal-to-noise near the edges of the spectrum. The high 
frequency component also varies slightly with spatial po- 
sition. The best correction is achieved if both object and 
atmospheric standard star have the same point spread func- 
tion (PSF), so that the pattern divides out during the flux 
calibration. The easiest way to ensure this with CGS4 was 
to smear the PSF slightly for both. The data from the sepa- 
rate beams were treated individually since the CVF pattern 
is seen to vary at the few % level between the two beams. 

The array was stepped six times across two pixels in 
order to fully sample the resolution element, and reduce the 
impact of bad pixels. These individual frames must be in- 
terleaved to give the final data and it is possible that any of 
sky transmission, atmospheric seeing or the position of the 
object in the narrow slit can change between frames. This 
is corrected for by looking for appropriate periodicity in the 
final spectra (ie a pattern repeating every 3 pixels) which 
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Figure 1. Spectra of the eight target stars and of the known Fell emitter IRC+10420 in the wavelength range showing HI Brll 
1.681/im and the 1.688/^m z 4 F9-c 4 Fg Fell line. Both the Fell line and HI Brll rest wavelengths are indicated by the dotted lines, and 
the previously unidentified feature we associate with Fel is also shown. The spectra are presented in the rest frame after correction for 
their published v^sr values for easier comparison. 
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Figure 1 - continued Spectra of the eight target stars and of the known Fell emitter IRC+10420 in the wavelength range showing HI 
Brl2 1.641//m and 1.644/im a 4 F 9 -a 4 D 7 [Fell]. Both the [Fell] line and HI Brl2 rest wavelengths are indicated by the dotted lines. The 
spectra are presented in the rest frame after correction for their published vlsr values for easier comparison. 
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can then be used to apply a correction factor to the data. 
Generally such corrections are small (<1%) but important 
in obtaining the best possible signal-to-noise ratio from the 
acquired data. 

Wavelength calibration was achieved using the telluric 
atmospheric absorption features present in all spectra. There 
are significantly more of these than either arc lines or night 
sky lines and hence they provide a much better map of the 
distortion of the wavelength scale and of the variation from 
object to object. We measured the wavelengths of the ab- 
sorption features from the much higher spectral resolution 
atmospheric absorption spectrum published by Hinkle, Wal- 
lace & Livingston (1995). Comparison of the different HI 
Brackett lines showed this procedure gave consistently good 
results, allowing us to compare the different lines with con- 
fidence. Our observations of NGC 7027 indicated our wave- 
length calibration should be good to within ±8kms~ . 

Finally, after objects are ratioed with the standard stars 
the two beams are combined to give the final flux calibrated 
spectrum. The final spectra are shown in Figure [1] 



3 RESULTS 

3.1 General Properties 

The hydrogen Brackett Br 11 and Br 12 lines were detected in 
all objects observed, as is evident in Figure[T]and Br7 is also 
present in all 5 sources we observed at that wavelength as 
shown in Figure[3] We also detected the 1.688/im z 4 Fg-c 4 Fg 
Fell in all sources, though in M8E it is very weak. Forbid- 
den [Fell] emission was also clearly present in BD+40°4124, 
M17SW IRS1, MWC 349A and S106 IRS. In S106 IRS the 
[Fell] emission is extended, with a strong velocity gradient 
along the full length of the slit. The bulk of the emission for 
our slit position angle comes from the north-eastern cone of 
the bipolar nebula (Solf & Carsenty 1982), and is blueshifted 
with regard to the systemic velocity. For each of these four 
objects the [Fell] line width was SOkms" 1 or less. The [Fell] 
line is a known tracer of fast shocks, such as seen in col- 
limated outflows near young stellar objects, and the line 
widths observed are perfectly consistent with this explana- 
tion. There is also a weak broad feature in both GL 490, 
and possibly VV Ser, that may be due to [Fell] but is offset 
to the blue from the systemic velocity (by about 60kms~ 
and lOkms -1 respectively). However the width of the line in 
these cases, about lOOkms -1 for VV Ser and 200kms _1 for 
GL 490, suggests a different origin from the other sources un- 
less we are seeing a one-sided jet coming towards us in these 
sources. These are the also two lowest luminosity sources 
we observed. We have been unable to identify any possi- 
ble alternative identification, but it seems probable that the 
line seen is not [Fell]. With these two exceptions, the for- 
bidden [Fell] emission has a very different line profile from 
the Fell line, and is always much narrower than any of the 
other lines observed. In addition, the [Fell] line is seen in 
NGC7027 whilst the fluorescent Fell is not, as expected for 
a purely nebular origin for the emission in this source. We 
can be confident therefore that the [Fell] line has an origin in 
the extended nebular gas around our sources when present, 
and Fell line has a different origin. We will not consider the 
[Fell] line further in this paper since we primarily observed 



it to demonstrate that the fluorescent Fell emission was not 
excited in a similar fashion to the [Fell] emission. 

An otherwise previously unidentified line at 
~1.68955/im is clearly visible in GL 490, M8E, MWC 297 
and MWC 349A, and possibly in S106 IRS. In S106 IRS 
the width and double peaked nature of the main Fell line 
blends with this feature, making the exact identification 
difficult, but there is excess red emission present that would 
naturally be explained if this unidentified line was also 
present here with a similar profile to that of Brl2 and 
Fell. Both Fel and CI lines are identified at 1.68959/zm in 
the solar spectrum (Ramsauer, Solanki & Biemont 1995). 
Observations of Fel 8387A in IRC+ 10420 (eg. Oudmaijer 
1998) show an almost identical line profile to the line shown 
here in that object (Figure Q]). We therefore prefer an 
identification with Fel for this line. Another unidentified 
line can be seen near 1.6385/im in IRC+10420 and possibly 
M8E, which again might be identified with a blend of two 
Fel lines (from the atomic data presented in Melendez 
and Barbuy 1999). This is certainly consistent with its 
appearance in only these two objects, but the profile is of 
poorer signal-to-noise in IRC+10420, which, together with 
the possible blending, makes it difficult to compare the 
profile of this line with that seen at 1.68959/^m. 

Figure [2] and Figure [3] present the key observational re- 
sults that we will interpret in the rest of this paper, namely 
a comparison of the profiles of the Fell, Brl2 and Bi'7 lines. 
The continuum has been subtracted from all the profiles. We 
note that weak Hel line emission is present in both S106 IRS 
and MWC 349A, on the blue wing of Brl2 and, more obvi- 
ously, Br7, but we have not attempted to correct for this. 
We have corrected the observed lines for extinction in order 
to make a detailed comparison. We adopt the extinction law 
from Cardelli, Clayton & Mathis (1989) to convert the liter- 
ature values in Table [1] to Ah and Ak, using the standard 
value of R — 3.1. The data shown in Figure [5] and Figure 
[3] have had these corrections applied. Finally we have nor- 
malised the plotted data such that the peak intensity of Brl2 
is given a value of unity, and all other corrected fluxes scaled 
relative to this. 

In order to compare the data with previously published 
work, and as a crude guide to the properties reflected by the 
observed line emission, we have also fitted all of the lines 
with a Gaussian. This was done to determine the systemic 
velocity as measured by that line, vlsr, and the approxi- 
mate full width at half maximum (FWHM). Where double 
peaked line profiles are present single Gaussian fits were used 
to measure line centre and FWHM, with the fit restricted 
to the line wings. This gives reasonably accurate values for 
vlsr but the line width is only accurate to approximately 
10%, and the quoted error reflects this. The full width at 
zero intensity (FWZI) is measured from the point where the 
line is 3a above the noise level, after smoothing the data 
by a factor of two, for direct comparison with the data pre- 
sented in Bunn et al. (1995). The results are shown in Table 
1 

Finally we note that the objects are split almost evenly 
between those where the overall general line shape is similar 
in Brl2 and Fell, and those which show significant differ- 
ences. However, the hydrogen lines are as broad or broader 
than the Fell line as measured at the half maxima points, 
and in BD+40°4124, GL 490, M17SW IRS1 and VV Ser the 
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Figure 2. Brl2 (solid) and Fell (dashed) line profiles after continuum subtraction. The velocity scale for all objects is in their own rest 
frame, after correction according to the nominal vlsr values given in Table 1. The spectra are normalised so that the peak flux of Brl2 
is set equal to 1, and the Fell line is scaled relative to this. 
© 0000 RAS, MNRAS 000, 000-000 
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Object Brl2 

vlsr FWHM FWZI Flux 
BD+40°4124 14± 1 320±30 ±510 1.61 ± 0.05 
GL 490 -24 ±1 192 ± 1 ±275 0.4 ± 0.02 

M17SW IRS1 14 ± 1 210 ± 20 -420/ + 330 0.3 ± 0.02 
M8E -37 ±1 116 ± 3 -600/ + 730 0.9 ± 0.05 

MWC 297 10 ±1 64 ± 1 -375/ + 400 12.4 ± 0.01 
MWC 349A 7±1 83 ± 8 ±330 24.0 ±0.1 
S106 IRS -6 ± 1 280 ± 30 ±420 0.39 ± 0.02 
VVSer -18 ± 1 380 ±2 ±420 0.52 ±0.1 



vlsr 
8±1 



FWHM 
233 ± 1 



S± 1 163 ± 1 



15 ± 1 97 ± 10 
-1 ± 1 185 ± 20 
-3 ± 1 320 ± 5 



Br7 
FWZI 
-500/ + 530 

±550 



Flux 

3.2 ±0.1 

1.3 ±0.01 



-300/ + 400 120 ± 1 

/ + 380 4.54 ±0.01 
-600/ ±700 2.1 ±0.1 



V LSR 

3±1 
-32 ± 1 
4± 1 
9±1 
6±1 
2±1 



Fell 
FWHM FWZI 
74 ± 1 ±100 



97 ± 1 
80 ± 1 
71 ± 2 
44 ± 1 
68 ±7 



19 ± 1 235 ± 20 
-6 ±1110 ±10 



Flux 
0.29 ±0.01 
0.1 ± 0.01 
0.1 ± 0.01 
0.02 ±0.01 
1.06 ±0.03 
13.50 ±0.05 
-200/ 0.13 ±0.02 
±80 0.07 ±0.01 



±80 
±80 
±55 
±50 
±100 



Table 2. Measured observed line parameters: all velocities are quoted in kms -1 , and all fluxes in units of 10~ 15 Wm _2 /im _1 . The FWZI 
missing for the blue wing of Br7 and the red wing of the Fell line in S106 IRS are due to strong residual line emission on the wings of 
the tabulated lines from other species. The errors quoted are the nominal errors on the fit, and do not include any error in the velocity 
calibration itself. 



difference is a factor of 2-3 as is easily seen in Figure [2] and 
Table [5] The hydrogen lines also have much higher veloc- 
ities in the line wings as measured at zero intensity, with 
the possible exception of S106 IRS. In all cases the veloc- 
ities observed are too large to be explained by gas excited 
in a standard HII region. There is a rough trend for more 
inclined objects to have broader Brl2 lines, as would be ex- 
pected when the motions in the disk plane are more along 
our line of sight (eg Home and Marsh 1986). 



3.2 Comparison of the hydrogen line profiles 

The comparison of velocity resolved hydrogen line ratios can 
potentially reveal the origin of the emission. If the gas is 
optically thin, this ratio will be a constant and will tend to- 
wards the case B ratio (eg Storey & Hummer 1995). Bunn 
et al. (1995) however showed it was possible to define the 
origin of the emission in the case where the lines are at least 
partially optically thick. Brll and Brl2 are too close to- 
gether to provide a useful comparison of line opacity (and 
as noted above we discount the Brll lines from analysis 
given the problems associated with its location near the 
edge of the observed spectra). We did observe BD±40°4124, 
M17SW IRS1, MWC 349A, S106 IRS and VV Ser in both 
Br7 and Brl2 lines. Figure [3] shows the Br7 and Brl2 lines 
for each of these objects, corrected for the inferred extinc- 
tion as described in Section 3.1, as well as the ratio of these 
lines. Away from the bright line cores, the main source of 
noise is actually non-Gaussian read-noise from the array. We 
have therefore truncated the line ratio plots at the velocity 
where the noise starts to dominate the ratio, but it is evi- 
dent from the over-plot of the lines themselves that the ratio 
remains roughly constant even beyond this. The presence of 
Hel emission is evident in MWC 349 and S106 IRS, at a sys- 
temic velocity of ~ — 200kms _1 relative to Br7, as a turn-up 
in the ratio at that point. 

The ratio in absorption cross-section between Br7 and 
Brl2 is 11.1, whilst that between Bra and Br7 is 15.9, and 
Bra and Pfy is 14.0. The comparison of Br7 and Brl2 is 
therefore a valuable test of the line formation mechanism in 
the same sense as the comparison of Bra and Pf7 or Br7 in 
the work of Bunn et al. (1995). For the two objects we have 
in common with Drew et al. (1993) and Bunn et al. (1995), 
the line ratios actually show very similar results near the 



systemic velocity (Figures 3(b) and 3(d) l. The earlier data 
have slightly lower spectral resolution and are also of lower 



signal-to-noise, which together can explain any differences 
seen. 

The optically thin case B value of the Br7/Brl2 ratio 
is ~ 5 ± 1 for values of the electron density in the range 
10 4 - 10 9 cm~ 3 and temperature in the range 3000 - 30000K 
(Storey & Hummer 1995). Simple theory indicates that for 
completely optically thick lines the ratio of the intensities 
of the two should scale simply as h/h = (A2/Ai) 4 (5i/5 , 2) 
in the Rayleigh- Jeans limit, where Si is the effective surface 
area of the emitting region of line 1. //the emitting region 
of the two lines has the same surface area then the limiting 
ratio is ~0.36 for Br7/Brl2. In practice we would expect the 
lower excitation Br7 line to emit from a larger volume, and 
hence larger surface area, so a value somewhat above this 
would be expected for gas which is close to optically thick. 
Figure [3] shows all five objects, with the exception of gas near 
the systemic velocity in S106 IRS, lies between the case B 
and optically thick limit, indicating opacity in at least one 
of the Brackett lines. Overall the trend in the ratios is very 
similar for four of the objects, namely a largely symmetrical 
profile that peaks near systemic. The exception is MWC 349, 
which has a very asymmetrical profile. We therefore defer a 
full discussion of a suitable model for MWC 349 to Section 

MM 

The generic features that any theory for the origin of 
the hydrogen line emission must be able to explain can be 
summarised as follows. First, we reiterate that all of our 
sources show all of the observed hydrogen lines in emis- 
sion, and the objects show evidence for nearly symmetric 
very broad line wings (Figures [5] and [3} • Secondly, many of 
the line ratios tend to be closer to the optically thick limit 
(Fsr-y I Fbt\2 ~ 1) in the line wings than in the core. Last, 
three of the objects show evidence for weak double peaked or 
self-absorbed line profiles, in either Brl2 or Br7 or both. In 
two of the sources the minimum is significantly offset from 
the systemic velocity: in BD±40°4124 it lies at ~ 21kms~ 1 , 
and in M17SW IRS1 at ~ 25kms~ 1 . In S106 IRS however 
it lies at ~ — Skms^ 1 . This however is within the error on 
our velocity scale, and hence is also in practice consistent 
with a small redshift. The first two are suggestive of self- 
absorption due to infalling material. The latter may still be 
self-absorption, or an intrinsic double peaked profile centred 
on the systemic velocity. These three objects are either at 
intermediate inclination angles or close to edge-on (Table[T|. 

These characteristics are shared by at least some lower 
mass pre-main sequence stars as well. Folha & Emerson 
(2001) present Br7 spectra of a large number of T Tauri 
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(c) MWC 349A (d) S106 IRS 

Figure 3. Br7 and Brl2 line profiles after continuum subtraction, and correction for extinction. The Brl2 line has been scaled as in Figure 
[2] and the Br7 lined scaled relative to this. Hel 7 3,1 G-4 3 ' 1 F emission is evident in MWC 349A between a velocity of ~ —250 — s> — lOOkrns -1 
on the Br7 spectrum. The same emission is also present in S106 IRS though only shows because of the asymmetry of the blue and red 
line wings to Br7. We have not attempted to correct the plotted ratio for this Hel emission. The ratio shown is for the line profiles after 
the corrections stated above have been applied. The ratio is smoothed by a factor of two relative to the line profiles, and truncated at 
the point where noise becomes dominant as noted in the text. The velocity scale is in the rest frame of the object as in Figure [2] 



stars. Symmetric broad wings and weak Br7 line self- 
absorption are common in the objects they class as type I 
under the scheme of Reipurth, Pedrosa & Lago (1996: type I 
objects in the original work refer to objects with symmetri- 
cal and unabsorbed Ha profiles). These are all objects that 
show greater evidence of veiling in the continuum, which is 
usually taken to indicate a high accretion rate. Objects with- 
out this veiling show less evidence for strong broad emission 
lines. It is worth noting that many of the strong T Tauri Br7 



line emitters still show P-Cygni style wind absorption, but 
only in the Hel lines that end in the metastable 2 1,3 5' levels 
(eg for the 1.083/mi 2 3 S - 2 3 P line, Edwards et al 2006). 
Many Herbig Ae/Be stars also show pronounced P-Cygni 
profiles for the Hel 1.083/mi line (Oudmaijer et al. 2011) as 
does MWC 297 (Drew et al. 1997), and another object in 
our sample, S106 IRS, has a P-Cygni profile in the 2.058/im 
Hel 2 1 S - 2 X P line (Drew et al. 1993). 

The natural explanation for the Hel profiles is an out- 
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flowing wind. The terminal wind velocity is then given in the 
usual fashion by the blueshifted absorption. For both of the 
objects in the sample under consideration here, S106 IRS 
(Drew et al. 1993) and MWC 297 (Drew et al. 1997), this 
gives a probable terminal velocity of ~ 200kms _1 . The wings 
on the hydrogen lines are still broader than this however (cf 
Figure [3]). In addition the wings have approximately the 
same full width at zero intensity in Brl2 as Br7. As noted 
by Bunn et al. (1995), if an accelerating wind is present the 
expectation is that the Br7/Brl2 line ratio should increase 
at high velocities: the faster gas lies physically farther out, in 
less dense gas, and is therefore optically thinner, and there- 
fore the full width at zero intensity should reduce for higher 
series members. None of our observed ratios match this. We 
can therefore rule out a standard hot star accelerating wind 
as the source of the emission. A decelerating (perhaps mass 
loaded) or constant velocity wind cannot be ruled out how- 
ever. 

The obvious alternative would be a variant on a disc 
model. However, this would be expected to show double 
peaked line emission if edge-on and relatively narrow emis- 
sion without broad wings if face-on for a pure disc model. 
A combination of wind and disc might appear to solve the 
problem, but such a geometry effectively masks much of the 
receding emission in stars that are near face-on and the sym- 
metry of the lines argues against this being the sole mech- 
anism responsible (see Section |3.4.5[) . An interesting theo- 
retical study with regard to the modelling of T Tauri stars, 
by Kurosawa, Romanova & Harries (2011), makes many of 
these same points about the difficulty of fitting profiles to the 
strong Hel absorption, Br7 emission sources. Their model 
includes potential emission from a stellar wind, a disc wind 
and magnetospheric accretion. They are able to reproduce 
the need for blueshifted absorption in Hel, and mostly sym- 
metric emission in the infrared hydrogen lines. However, the 
lines are too narrow, and at high mass accretion rate the hy- 
drogen lines start to show inverse P-Cygni absorption from 



the infalling gas, suggesting some extra component is re- 
quired to fully explain the broad wings in the hydrogen lines. 
Notably, the disc wind also must have a large opening an- 
gle to avoid significant self absorption even at relatively low 
accretion rates, and the same is likely to be true for our 
sources. 

This raises the possibility that additional line broaden- 
ing may be present in the hydrogen lines. Typical densities 
of the inner regions of an outflowing wind in these sources 
are naively of order > 10 12 cm -3 assuming the mass loss rate 
is > 10 -6 Mq /yr, the observed wind velocity of ~ 200kms _1 , 
and typical main sequence OB star radii. The radius may ac- 
tually be an underestimate since models predict the stars are 
swollen until they reach a final main sequence configuration 
(eg Hosokawa & Omukai 2009) , consistent with the low stel- 
lar wind speeds we see. However, the mass loss rates are also 
likely to be substantially higher, so the density estimate is 
relatively unchanged. This density is large enough that Stark 
broadening can become significant for some lines (cf Muze- 
rolle et al 2001 for the case of T Tauri stars). Repolust et al. 
(2005) modelled the Stark broadened Brackett line profiles 
for main sequence OB stars. Although the winds present in 
main sequence stars are very different, the density is sim- 
ilar to the estimate given above for our simple outflowing 
wind, and the principle of whether broadening occurs under 
such circumstances is identical. Figure 3 in Repolust et al. 
shows that in this density range wings are easily generated 
by Stark broadening for BrlO, with broad wings similar to 
that seen in our data. Electron scattering can also be appar- 
ent with such high densities. With T e ~ 10000K, line wings 
comparable to the thermal velocity width of the electrons 
of 550kms _1 is possible. We know scattering is significant in 
some Herbig Be stars, including MWC 297 (eg Oudmaijer & 
Drew 1999). Another key advantage of a broadening mech- 
anism is that it can help to fill in any depression in the line 
caused by line self-absorption due to shadowing from a disk 
(eg see the schematic in Malbet et al. 2007 with regard to 
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MWC 297). Finally, broadening decouples the line ratio at 
high velocity from the need for an optically thick emission 
region at that velocity. 

In principle we can even tell which broadening mecha- 
nism is actually present. A simple comparison of the ratio of 
the intensity of the line wing to the peak intensity will show 
stronger wings for higher series transition if Stark broaden- 
ing is present, since these are the states most affected by per- 
turbations due to neighbouring ions. We would expect the 
opposite to be true for electron scattering, where essentially 
it is the population of the excited state that determines the 
likelihood of scattering. Thus a line like Ha is more likely to 
have wings enhanced by scattering, and the relatively high 
Br series lines we are studying are more susceptible to Stark 
broadening. Stark broadening should also tend, in the high 
velocity limit, to wings that scale as v~ 5 ^ 2 (Repolust et al. 
2005). The core of the line, as shown clearly in the data of 
Repolust et al., should also be almost flat in a log Fa ver- 
sus logn plot. Scattering has been modelled previously in 
the young stellar object Lk Ha 101 by Hamann and Pers- 
son (1989), using the prescription of Castor, Smith and van 
Blerkom (1970). In this model the scattering all occurs in a 
slab that lies just above the line source (the stellar wind), 
and is optically thin to Thomson scattering. Laor (2006) de- 
scribes a model for scattering in the broad line region in ac- 
tive galaxies, with a similar requirement on the opacity but 
an isotropic illumination of the electrons, rather than a one 
sided illumination. The resultant functional forms for the 
scattered flux are different, but both give profiles that are 
roughly exponential for likely values of the two key scatter- 
ing parameters (opacity and electron velocity). In principle, 
the Castor et al. model can be used to measure these scat- 
tering parameters, and hence also quantities such as density 
at the base of the wind. In practice, this has not proved pos- 
sible as we shall discuss in our analysis of the best fit to the 
line profiles in Section 3.4, and in Section 4. 

Finally, as noted by Drew et al. (1993), an optically thin 
nebular component may be required in some cases in addi- 
tion, of which S106 IRS is a clear example given it powers 
an extended bipolar HII region (all three Brackett lines pre- 
sented are visible across an extended region on the long slit 
spectrum, with lines of width ~ 30kms _1 ). In that instance 
the line ratio at systemic velocity is actually above the case 
B value. This most likely indicates that the optically thin 
nebular gas actually has lower extinction than S106 IRS 
itself (ie we have over-corrected for the extinction in this 
component), which is consistent with the fact the nebula is 
visible in the optical and the star itself is extremely faint. 

Overall therefore, broadening gives a natural explana- 
tion for the wings in all our sources. If we accept that factor, 
then a decelerating, constant wind or large angle disc wind 
are viable explanations for the core of the lines. This is prob- 
ably a natural explanation for all massive and intermediate 
mass young stellar objects since none of our objects are se- 
lected to be in any way "special" . 

3.3 Modelling of the Fell line with a disc 

The fluorescent Fell emission is a clear tracer of emission 
from a disk in classical Be stars (Carciofi & Bjorkman 2006). 
Detailed models have been derived for observations of the 
optical Fell lines in Be stars (eg Arias et al 2007), show- 



ing that even quite complex line profiles can be explained in 
terms of emission from a relatively narrow ring. It therefore 
seems natural to consider whether a disc model can explain 
these relatively younger objects. Currently, the 1.688/im Fell 
line is not yet included in any radiative transfer model. Nor 
do we understand its source function well enough, given its 
potential origin as a laser (Johansson & Letokhov 2007), to 
reliably model even in the fashion of Arias et al. We there- 
fore rule out detailed calculations at this point, and leave 
that to future work. Instead, we adopt a "proof of concept 
approach" and attempt to fit the observed line profiles with 
a simple accretion disc model but without a detailed source 
function for the line emission, or attempting to enforce a 
correct temperature structure for the emission zone. Specif- 
ically, we use the constraints given by other observations to 
fix the mass of the star and the inclination (as quoted in 
Table [T]) , and then test how well a disc can fit the observed 
line profile. 

The model in question consists of a Keplerian disc 
within which the line flux per unit area is proportional to 
the local surface density. The structure of the rotating disc 
is given by the standard Q-disc model - as described by, eg, 
Pringle (1981). The accretion rate and value of a are not 
fitted in this process, since the model is being fitted to a 
normalised line profile, and thus the actual surface density 
is unimportant. Furthermore, the variation of surface den- 
sity with radius is insensitive to both a and the accretion 
rate. As a result, the only free parameters are the intrinsic 
line width and the outer radius of the emitting area. The 
latter is poorly constrained for most of the objects however 
since they exhibit single peaked emission. In that instance 
the outer radius is a lower limit set by the model line width, 
since we clearly cannot constrain this parameter strongly in 
the case of a near face-on inclination. 

Turbulent line broadening is added to the disc model 
in order to achieve the final fits. This is represented by a 
Gaussian, and the width is given in Table No other line 
broadening mechanism is considered as we do not expect any 
to occur - even thermal broadening is likely to be insignifi- 
cant for iron. Indeed, as noted by Johansson and Letokhov 
(2007), the conditions under which this Fell line is emitted 
are likely to lead to mild lasing, which actually narrows the 
natural line profile. We also do not expect the opacity in 
this transition to be significant. 

The model is used to determine the line profiles for a 
given set of input parameters. The model profiles are then 
smoothed to the instrumental resolution, before being com- 
pared to the observed profiles. To find the best fitting pro- 
files, the outer radius and line width are allowed to vary be- 
tween 1-100AU and l-75kms~ respectively and the Powell 
method (as implemented in idl) is used to minimize the 
residual of the fit. We present a comparison of the best fit- 
ting model profiles and the data in Figure [4] A detailed 
case-by-case discussion of the best-fitting parameters is de- 
ferred to the next section, but it is clear that, in general, 
the observed Fell 1.688/^m line profiles can be modelled by 
emission from circumstellar discs of an appropriate size. The 
major exception to this is S106 IRS, which we discuss in de- 
tail in Section [3A7l Figure [2] shows that GL 490, MWC 297 
and MWC 349A also appear similar in Fell and Brl2 if we 
ignore the broad Brl2 line wings. A disc model can also 
partly explain the Brl2 emission in all of these sources as 
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Figure 4. Fell emission line profiles of each object and the associated best-fitting disc model profiles (dashed line). In VV Ser the "spike" 
evident in Figure [2 (h)| has been clipped to allow a better disc fit to the underlying line profile. 
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6.1 


3.56 
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21.6 
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4.1 


GL 490 


8.6 


4.28 


30 


5.2 


75.0 
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20 


15.6 


44.6 
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15 


3.7 
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67 
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68.8 


0.7 
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10.0 


6.12 


15 


10.5 


36.5 


0.5 


MWC 349A 


26.0 


8.54 


85 


31.0 


26.4 


11 
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21.7 


7.51 


80 


1.8 


26.8 


1.0 


VV Ser 


3.8 


2.81 


85 


2.0 


50.1 


0.8 



Table 3. Parameters of the best fitting models. The stellar radii and masses are taken from the compilation of Harmanec (1988). For 
M17SW IRS1 we show two fits to the data, one reflecting the inclination inferred by Follert et al. (2010), the other for a more face-on 
orientation. The fit shown in Figure [4(c) | is the former, but there is essentially little difference between the two. The other inclinations 
are adopted from Table [l] 



discussed below, consistent with the results from fitting the 
iron line. A disc model does not fit the HI lines from the 
other sources well at all. 



3.4 Comments on individual objects 

3.4.1 BD+40°4124 

The Ha line in this source shows an asymmetric double 
peak (Vink et al. 2002), closer to line self-absorption than 
to evidence for a rotating disk. The velocity separation be- 
tween the peaks is ~120kms _1 , compared to only about 
~ 40kms _1 as seen in our Br7 spectrum. This is also consis- 
tent with the self-absorption picture since we would expect 
Ha to be more optically thick at velocities near, and to the 
blue of, systemic compared to Brl2 if a dense wind is present 
(as is evidenced by Figure 3(a) I. Eisner et al. (2003) find an 



inner radius of 0.73au from modelling their continuum inter- 
ferometric data with a flat dusty disc. Any disc within this 
radius must be dust-free. Vink et al. (2002) find that their 
Ha polarisation data can be explained by a straightforward 
line depolarisation effect, where the continuum light is scat- 
tered by a small scale disc, whereas the emission line itself 
arises on scales larger than the disc, likely in a wind. The 
small scale disc could be the extension of the dusty disc de- 
tected by Eisner et al. down to the star. The line wing/peak 
ratio is roughly the same for Br7 and Ha, but larger for 



Brl2 (cf Figure 3(a) I, suggestive of Stark broadening. The 



line wings can be fitted equally well by both a power-law, 
with exponent ~ —2.0 ± 0.1, and an exponential (Figure 
[5|. We also explicitly fitted the Castor et al (1970) model. 
A formally good fit is achieved, but the key parameters in 
this case, the width of the line core, the electron tempera- 
ture (width) and opacity, are highly degenerate. Effectively 
there is little difference between a very broad intrinsic line 
and a highly scattered one. The wing-peak ratio for the Brl2 
and Br7 data, and the closeness of the power-law exponent 
to the canonical —2.5, both suggest that Stark broadening 
is present however. The Fell emission in this object is not 
fitted perfectly by a disc model, but the profile is clearly 
very different from the Brackett lines. In this instance it is 
plausible that the hydrogen lines arise in a large scale wind, 
whereas the Fell is mostly tracing the inner disk. 



3.4.2 GL 490 

Schreyer et al. (2006) detect a relatively large-scale (~ 
1500au) molecular disc in approximate Keplerian rotation 
from which they derive an inclination of ~ 30°. The posi- 
tion angle for this disc is the same as that found at 2cm by 
Campbell, Persson & McGregor (1986), and both lie orthog- 
onal to the direction of the observed large scale molecular 
outflow observed by Mitchell et al. (1995), similar to other 
equatorial wind systems (eg Hoare 2006). Oudmaijer, Drew 
& Vink (2005) also inferred the presence of a small scale disc 
from near infrared spectropolarimetry. Oudmaijer et al. did 
however require the initial line emitting material to lie inter- 
nal to the scattering disc, and that the disc did not extend 
down to the stellar surface. 

Our Brl2 profile looks very similar to the Br7 profile 
in Bunn et al. (1995). Bunn et al. also infer a relatively flat 
Bra/Br7 ratio in GL 490 out to velocities ~ ±100kms -1 , be- 
fore an upturn is seen. However, the evidence for a turn-up 
in these earlier CGS4 data may be compromised by poorer 
continuum subtraction, since the data were taken with a 
smaller array and at lower spectral resolution, spatial reso- 
lution and free spectral range. Comparison of our Brl2 pro- 
file with the Br7 and Bra profiles in Bunn et al. suggests 
that the broad line wings are relatively stronger in Brl2 
which slightly favours Stark broadening. There is a good fit 
to the Castor et al model for Br 12 as well, and in this in- 
stance the parameters are well defined with r ~ 0.3 ± 0.05 
and T e ~ 340 ± 30K. The very low value of the electron 
temperature implied by the fit is a concern, but we do not 
have sufficient signal-to-noise to determine if the wings ex- 
tend further in velocity. The line wings do not follow the 
expected power law for Stark broadening either, having an 
exponent of only ~ — 1.2±0.2. In this instance, the relatively 
poor signal-to-noise prevents a truly firm conclusion. 

The Brackett lines differ from the Fell line largely in 
a strong red asymmetry, and stronger evidence for broad 
line wings in the hydrogen lines. The Fell emission can be 
successfully fitted with a Keplerian disc. The inferred lower 
limit to the outer radius is much smaller than the size of the 
structure seen in the radio map which is several hundred AU 
across. Given the similarity of the line core in this object 
between Brl2 and Fell it is possible that a disc may be a 
suitable explanation for the hydrogen line emission as well, 
especially given the orientation of the radio emission. There 
is no self-absorption in any of the lines observed either here 
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Velocity (km/s) Velocity (km/s) 

(a) Power-law fits (b) Exponential fits 

Figure 5. Fits to the Brl2 line for either a power law or an exponential. The fits are shown as the thick dashed (red) line. The data have 
been binned into lOkms -1 bins, and folded around the systemic velocity. The fits are made only to the region between 100 and 400kms — 1 . 
The objects shown in each plot are (from top down): VV Ser, S106 IRS, MWC 297, M8E, M17SW IRS, GL 490 and BD+40°4124. 



or in Bunn et al. so there is no requirement for a strict wind 
component to be present. 



3.4.3 M17SW IRS1 

M17SW IRS1 is a young binary star located just beyond 
the edge of the M17 HII region and is also known as the 
Kleinmann- Wright object (Kleinmann & Wright 1973). Es- 
timates of its luminosity range from about 5000 L© (Chini 
et al. 2004) to 50000 L (Follert et al. 2010). The latter is 
more consistent with the presence of Hel 2.058^tm emission 
as reported by both Porter et al. (1998) and Hanson et al. 
(1997). 

The deredenned flux ratio is very similar to that pre- 
sented by Bunn, Drew & Hoare, and as noted by them 
is rather similar to that for BD+40°4124. Br7 does show 
weak red asymmetric self-absorption which might suggest 
the presence of infalling ionised gas. 

The Brl2 and Fell emission have very different profiles, 
again in a similar fashion to BD+40°4124. The Fell data can 
be explained by a disc, at both a near edge-on inclination as 
measured by Follert et al. (2010), as well as nearer to face-on 
(cf Table [3]). However, an edge-on disc requires a very large 
spread in the Fell emission region, almost out to the full 
extent of the mid infrared emission region seen by Follert 
et al. Without further information however we cannot rule 
out either possibility. The degeneracy between inclination 
and outer radius of disc is of course a well known problem 
in this situation, but M17SW IRS1 is the only source in 
which the result of fixing the inclination at the value given 
in Table [T] leads to results that are then difficult to interpret 
in light of other data. The most likely model for this source 
is the same as discussed for BD+40°4124, namely a disc 
giving rise to the Fell emission and a wind the Brackett 
line emission. The wings on Brl2 are relatively weak in this 
source, but the power law profile of the wing is ~ —1.8 ± 
0.2, suggesting that Stark broadening may be present. A 
formally good fit to the Castor et al model can be achieved 
but the values of both opacity and electron temperature are 
again highly degenerate with the core line width, and no 
meaningful physical parameters can be obtained. 



3.4.4 M8E 

The Brl2 line in M8E is best described as a relatively narrow 
Gaussian superposed on a much broader "pedestal". This 
is qualitatively different from the Bra and Br7 profiles in 
Bunn, Drew & Hoare. The line core is similar to Br7, but 
the higher velocity line wings are less evident. Bra has broad 
wings but the profiles are not a match. It is possible the 
greater wavelength coverage of the data presented here may 
partially explain this difference. It cannot be explained easily 
with simple line broadening mechanisms outlined in Section 
13.21 since there is no simple trend with quantum number n. 
An exponential fit is a better match to the broad wings than 
a power-law, and the latter also has a very low exponent of 
— 1.1±0.1. Neither matches at the highest velocities however. 
The Castor et al model is formally a very poor fit. 

One other notable aspect of our data is that the line 
centre velocity of the Brackett lines are significantly off- 
set with respect to systemic. As shown by Table [2] these 
lines have vlsr ~ — 37kms _1 , whereas the Fell line has 
vlsr ~ 9kms _1 , which is in agreement with the measured 
molecular gas value, and the Bra and Br7 data of Bunn et 
al. (1995). The fact that both Brll and Brl2 lines agree 
in velocity rules out a simple wavelength error in part of 
the data since they are observed in different echelle settings. 
Bunn et al. infer a relatively fiat Bra/Br7 optically thick ra- 
tio in M8E out to velocities ~ ±100kms _1 , with only weak, 
if any, evidence of a turn-up at higher velocity. 

The Fell line can be modelled by a disc, but the rela- 
tively weak line emission is noisy and hence does not partic- 
ularly constrain other models. Linz et al. (2009) can only set 
an upper limit on the size of the disc in the mid-infrared of 
~100AU, which is consistent with our Fell disc model. The 
similarity of the core of the Brl2 line and the Fell line, if 
we ignore the velocity offset mentioned above, suggests both 
have a common origin. Overall though the difference in ve- 
locity, the unusual "pedestal" on Brl2 and the fact that the 
line widths do not simply scale with n all make it difficult 
to ascribe an origin to the emission. 
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Figure 6. The best fitting Castor et al. (1970) model to MWC297. 



3.4.5 MWC 297 

The inclination of this source is the subject of debate. Drew 
et al. (1997) find a near edge-on inclination is necessary to 
explain the observed optical absorption line data, since oth- 
erwise the star is near break-up velocity. By contrast all of 
the interferometry data is best modelled by a near face-on 
orientation instead. An edge-on model is also more likely 
to lead to double-peaked line profiles if an ionised disc is a 
contributor. For consistency with the other objects studied 
here, where we have adopted inclinations from interferomet- 
ric data, we do the same here. 

The Bra/Br7 ratio presented in Murdoch & Drew 
(1994) is similar to the Br7/Brl2 ratio of BD+40°4124 or 
M17SW IRS1 shown in Figure [3] The line widths we mea- 
sure are similar to the Bra width as reported in Drew et 
al. (1997). The FWZI of Brl2 is closer to that seen in Ha, 
though, as is clear from Drew et al., the latter has a much 
larger FWHM. The Fell profile does not show the same 
broad wings. Most published spectra of MWC 297 show 
symmetric Brackett line profiles. The exception is Malbet 
et al. (2007), where a clear blue-shifted asymmetric self- 
absorption is present in Br7. The higher spectral resolution 
data presented by Murdoch & Drew (1994) or Weigelt et al. 
(2011) do not show this, which suggests variability is present. 
There is a clear P-Cygni effect seen in the 1.083/im 2 3 S-2 3 P 
Hel line (Drew et al. 1997), suggesting a wind is present, 
which for a relatively face-on orientation must either be a 
polar wind or a disc-wind with a large opening angle. 

Drew et al. present a 5GHz radio map showing resolved 
radio emission from this source on the scale of ~ 50/1(7, and 
both the Malbet et al. and Weigelt et al. interferometry also 
show the core Br7 emission is more extended than the K- 
band continuum light. The higher spectral resolution data 
of Weigelt et al. show the visibility is essentially the same as 
the continuum beyond about 60kms _1 from the line centre. 
The Brl2 wings is better fitted by an exponential in this 
case. Notably this is the only object in which the line core 
clearly has greater flux than a simple extrapolation of the 



exponential fit might suggest (Figure 5(b) I, as we would ex- 



pect if the exponential wings were light scattered from that 
core. A power-law fit by comparison again has a rather low 
exponent of —1.3 ± 0.1, and the flat inner line core from 
Stark broadening seen in the Repolust et al. data, and ob- 
jects such as BD+40°4124, is absent here. Comparison of 
the strength of the line wings to peak intensity in the hy- 
drogen recombination lines in the data presented here and 
in Drew et al. (1997) indicates that the wings grow stronger 
for the lower n transitions in a fairly uniform way. This is 
strongly suggestive of electron scattering. The fact that the 
wings and continuum appear to have the same visibility in 
the Weigelt et al. data suggests this scattering occurs on 
small scales. 

We consider the scattering model in more detail for 
MWC 297 since this source shows the best evidence that it is 
present. The best fit of the Castor et al model has a reduced 
X 2 = 2, indicative of the fact that the line is slightly asym- 
metric. We show the fit for this source in Figure [6] In this 
case the opacity to Thomson scattering is greater than unity, 
so the conditions on the actual model are not met. The im- 
plied electron temperature is also very low (T e ~ 170± 10K) 
and not physically meaningful. The match of the Castor et 
al profile in a case where the intrinsic line core itself is very 
narrow does argue for scattering being the dominant mech- 
anism, but that this model itself is insufficient to accurately 
estimate key physical parameters related to the scattering. 

If we accept a mostly face-on inclination for this source 
then our Fell profile can be fitted well by a disc model, with 
an outer radius that is well within the extent of the ionised 
gas as seen in the radio map. Acke et al. (2008) struggled 
to fit their combined near and mid-infrared interferometry 
with a pure accretion disc model, suggesting that more than 
one component was present. Malbet et al. suggest a combi- 
nation of disc and a weak outflowing wind can explain their 
near infrared interferometry and the observed P-Cygni pro- 
file they find in Br7, but they struggle to find a good fit 
to the red wing of the Br7 line because their model pre- 
dicts disc shadowing should be present where none is seen. 
The same is true of the models applied to the high spec- 
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tral resolution observations of Weigelt et al. (2011) as well 
(which otherwise reach a very similar conclusion to Malbet 
et al). However the relatively clear evidence for significant 
scattering in the line wings largely reduces this problem. In 
a face-on orientation, the lack of a clear line effect in the 
spectropolarimetry of Oudmaijer & Drew (1999) is not sig- 
nificant, since an asymmetry in the source structure as seen 
on the sky is required to give rise to such a signature, and 
this is naturally not present. A counter argument to this is 
that Weigelt et al. find evidence from the differential phase 
for a rotational signature in Br7, which of course we would 
not see in an almost completely face-on situation. An inter- 
mediate inclination, as suggested by Acke et al. (2008) as 
a compromise, may help to reconcile all of these pieces of 
evidence. 

3.4.6 MWC 349A 

MWC 349A clearly shows different line profiles and ratios 
to our other sources. The precise evolutionary status of 
MWC 349A is not known (see, eg, Hofmann et al. 2002 for a 
discussion), nor is an exact spectral type known. The main 
argument against it being young is that it lacks significant 
far-infrared thermal dust emission, with the dust emission 
peaking in the mid-infrared instead. The strong emission 
line spectrum and bolometric luminosity are also consistent 
with an identification as an evolved B[e] supergiant, though 
it also lies at the upper end of properties of known Herbig Be 
stars. The luminosity is 6 x 10 4 L0 (adapted from Cohen et 
al. 1985 using the distance from Meyer et al. 2002). This is 
consistent with the presence of Hel lines in the spectrum (eg 
Hamann & Simon 1986). MWC 349A shows strong hydro- 
gen maser recombination lines in the mm and radio regime 
(eg Weintroub et al. 2008), consistent with the edge-on in- 
clination found from interferometry. 

Fell emission from similar energy levels to the 1.688/im 
line also show double peaked profiles (eg the 2.089/xm Fell 
line in Hamann & Simon 1986, the optical lines in Hamann 
& Simon 1988) . These lines have been interpreted as having 
an origin in a disc (Hamann & Simon 1988). Higher excita- 
tion Fell lines appear more like the Br7 profile we see here, 
with a distinct line asymmetry (Hamann & Simon 1988). 
MWC 349 A is the object where our data have the highest 
signal-to-noise, and therefore small departures from a Ke- 
plerian disc are easier to identify. This is clear in our Fell 
disc fit, which formally is very poor. However, this is mostly 
due to a weak blue-red asymmetry in the observed data, 
which cannot be modelled within the context of our very 
simple prescription. Kraus et al. (2000) also had some dif- 
ficulty in fitting the v=2-0 CO bandhead with a Keplerian 
disc profile. A disc is clearly present in this source, however, 
given the maser emission. It may simply be that the disc 
is not completely Keplerian as also suggested by the maser 
observations of Weintroub et al. (2008). Danchi et al. (2001) 
derived an outer radius of about 50AU from near-infrared 
interferometry, consistent with what we find from modelling 
the Fell data, and the origin of the mm/sub-mm masers 
(Thum et al. 1994), though all are rather larger discs than 
the best fit that Kraus et al. find for their disc model. 

The HI Br7 profile also appears very similar to that 
shown in Hamann and Simon (1986), and shows what is 
likely to be weak blue self-absorption or a very asymmetric 



disk. The broad wings can again be ascribed to electron scat- 
tering (Meyer et al. 2002), in agreement with the spectropo- 
larimetry of Oudmaijer et al. (2005). The hydrogen lines 
have generally been interpreted as the combination of emis- 
sion from a slow stellar wind and a rotating disc (eg Hamann 
& Simon, 1986, 1988). This wind must be on larger scales 
than the disc modelled by Danchi et al., since otherwise the 
velocities would be anomalously low. Another possibility is 
a disc-wind since the low expansion velocity evident in the 
hydrogen lines naturally arises in a model where the wind 
is launched from the disc at some distance from the star (eg 
Gordon 2003, Hollenbach et al 1994, Sim et al. 2005). 

The Br7/Brl2 line ratio plot shows a distinctive pattern 
different from any of either our sources or those presented 
in Bunn et al. It is however consistent with self-absorption 
in the lines due to an outflowing wind in the sense that 
the opacity increases uniformly towards the blue. We also 
considered the ratio after the individual line profiles had 
been cross-correlated, in order to compensate for any resid- 
ual velocity calibration errors that might be present, which 
could have an effect here given the steepness of the pro- 
file. The result shows a slightly flatter line ratio, but still 
with a distinct slope in the same sense as shown here. Kraus 
et al. found that the high series Pfund lines were singly 
peaked with widths ~ 50kms _1 . This suggests these lines 
arise in optically thin gas far from the star in the disc-wind, 
or in the nebula. However such a large scale optically thin 
Pfund component should be reflected in both the Balmer 
and Brackett lines as well, and is not. The mid infrared in- 
terferometry presented by Quirrenbach, Albrecht & Tubbs 
(2006) also shows the high excitation mid infrared hydrogen 
lines are less extended than the dust continuum emission. 
Instead the Pfund and higher series lines must, on energetic 
grounds, arise near the central source, and the single peaked 
nature implies optically thin, approximately spherical, wind 
emission. It is probable therefore that the infrared recom- 
bination lines largely arise from a combination of disc-wind 
and spherical wind. 

3.4.7 S106 IRS 

S106 IRS is the exciting star of the Sh 2-106 HII region (see, 
eg, Saito et al. 2009). The HII region is best modelled as a 
bipolar flow which is almost oriented in the plane of the sky 
(eg Solf & Carsenty 1982). S106 IRS has a variety of reported 
extinction values in the literature. The value given in Table[T] 
is from a fit to the silicate absorption spectrum as measured 
by ISO in a large beam (van den Ancker et al. 2000) which 
will be dominated by the emission from S106 IRS. The fits 
to the CO bandhead by Chandler et al. (1995) shows that a 
disc of radius 1-2AU is likely to be present. 

The Br7 and Fell lines in S106 IRS are remarkably sim- 
ilar (Figure 2(g) I, though Br7 has a larger FWZI than the 



Fell line. Perhaps the greatest surprise is that despite the 
near edge-on nature of this source, a disc model is a very 
poor fit (Figure ID). We observed the Fell line in S106 IRS 
at two epochs separated by 378 days. The data shown in 



Figure 2(g) are the average of these observations. There was 



no evidence that the HI Brll line varied between these ob- 
servations but there is some tentative evidence of a change 
in the Fell line. This is shown in Figure [7] where we plot the 
data from the two epochs, after slight smoothing to improve 
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the signal-to-noise. There is a clear dip in the data at the 
velocity of maximum absorption for the second epoch. We 
can understand both the line profile and the variability if the 
Fell line in this case does arise in a clumpy (hence variable) 
wind. The remaining differences in the line profile between 
the Brf2 and Fell line then are due largely to the differing 
opacities in the two lines. We already know that Sf06 IRS 
shows strong P-Cygni absorption in the Hel 2 1 S-2 1 P tran- 
sition at 2.058/^m (Drew et al. f993), indicative of a stellar 
wind with a high mass loss rate (see also Felli et al. f984). 

The Br7/Brf2 ratio plot for this source is similar to the 
Bra/Br7 ratio presented in Drew et al. (f 993) but with bet- 
ter signal-to-noise at high velocity. Drew et al. found an up- 
turn at vlsr ~ — 200kms~ 1 in their Bra/Br7 ratio, which 
they interpret as evidence that Br7 is becoming optically 
thin. However we see no such evidence. Drew et al. tried 
to deblend the prominent 5 1,3 G-4 1,2 D Hel line from Bra at 
around this wavelength, whereas we have not done the same 
for the 2.16475/zm 7 3,1 G-4 3,1 F Hel lines which may explain 
the difference, though it should in theory enhance the up- 
turn in our data and suppress it in theirs. We therefore ac- 
cept our line ratio as indicative for the source as a whole. A 
broadening mechanism is therefore required to explain the 
line wings as in other sources. The Brf 2 wings have a power 
law profile with exponent —2.0 ± 0.2, similar to the case of 
BD+40°4f24. An exponential or the Castor et al model give 
a slightly poorer fit. Since our Castor et al model does not 
match the central absorption we do not consider the fitting 
parameters further for this source. The enhanced wings seen 
in the Brackett lines appear more prominently in Brl2 than 
Bra however, more suggestive of Stark broadening. 

An almost unresolved nebular core further complicates 
the picture for our Br7/Brf2 ratio. The nebular line arises 
from a much larger volume than the rest of the emission, and 
should have a lower extinction correction applied since the 
HII region is optically visible. This explains why the ratio 



plotted in Figure 3(d) exceeds the nominal case B value. 

Overall most of the lines in this source can be explained 
by a dense wind (as in Drew et al. 2003) which must be 
primarily along the equatorial direction to avoid veiling the 
nebula from the star. In this instance, the Fell emission may 
also arise from the disc-wind, given the asymmetric wind- 
like profile it has. 

34.8 VV Ser 

Blake & Boogert (2004) observed the M band fundamental 
CO transition in this source. Their results show a clear dou- 
ble peaked profile which they model as a disc with i = 85°, 
consistent with the interferometry of Eisner et al. (2003). 
Notably the blue peak of their CO profile is enhanced with a 
narrow emission component in the same fashion as our Fell 
profile. This is suggestive of a localised spot of enhanced 
emission, which must have an origin in a source with con- 
stant velocity (ie not the disc), since our data and the Blake 
& Boogert data are separated by 2-3 years. We do not quite 
resolve the blue and red peak in this source despite the al- 
most edge-on nature. The disc model fits well despite this 
if we ignore the extra "spike". Notably, Eisner et al. (2003) 
suggest that a disc of similarly small size is present in their 
interferometry, though they favour one with a puffed up in- 
ner geometry. Pontopiddan, Blake & Smette (2011) model 



the CO as arising at 3.5au, outside the emission region for 
the Tf-band continuum. VV Ser is observed to have a disk 
shadow, consistent with this modelling (see discussion in 
Pontopiddan et al. 2011). 

The Ha profile is asymmetric double peaked (Reipurth, 
Pedrosa & Lago, 1996) with a velocity separation of 
250kms _1 between the red and blue peaks, again sugges- 
tive of a weak P-Cygni profile. Our data lack such obvious 
P-Cygni effects, but the ratio shown in Figure 13.11 is oth- 
erwise similar to most of the other objects. An exponential 
fit to the line wings is marginally preferred. A power-law fit 
gives an exponent that is again rather shallow at — 1.5±0.1. 
The line wings are stronger relative to the peak in the lower 
series lines for VV Ser (eg Figure 3(e)), which is suggestive 
of a scattering origin as well. The Castor et al model gives a 
formally good fit but the opacity and electron temperature 
are effectively undefined in the result as was the case for 
other sources due to degeneracy in the best fitting model. 
A larger wavelength coverage would perhaps have helped 
identify the source of the broad line wings in this source. 
A wind origin for the remaining hydrogen emission seems 
likely given the very different Brl2 and Fell profiles, with 



the latter almost certainly coming from a disc (Figure 2(h) I 



4 CONCLUSIONS 

Our results show clearly that a wide variety of physical en- 
vironments can be traced using near infrared spectral lines. 
It is clear that the combination of HI, Fell and, where avail- 
able, Hel lines can help to disentangle the otherwise complex 
line profiles observed. Hydrogen line ratios are particularly 
powerful in showing the need for a line broadening mech- 
anism in all objects expect MWC 349. Specifically, data 
with sufficient spectral resolution and wavelength coverage 
allowed us to identify a likely origin in Stark broadening in 
the wings of hydrogen lines in BD+40°4124, M17SW IRS 
and S106 IRS, and electron scattered wings in MWC 297, 
and probably VV Ser. Better signal-to-noise, as available 
with larger telescopes, and an even larger free spectral range 
would have allowed us to be more conclusive about the origin 
of the wings in both GL490 and M8E as well. 

One important consequence of the presence of strong 
line broadening mechanisms in the high series infrared hy- 
drogen lines is that, on their own, they are poor tracers of 
a wind or a disc. A better tracer of a wind are undoubt- 
edly the Hel lines from metastable levels, although, unfor- 
tunately, not every young stellar object shows this emission. 
The hydrogen lines do however provide a means to test for 
the presence of very dense circumstellar material (through 
broadening). 

We also attempted to utilise the model of Castor et 
al. (1970) to estimate physical scattering parameters, as- 
suming that was the dominant mechanism. The fits were 
either highly degenerate in the parameters where the core 
line width is close to matching that required for the wings 
(eg BD+40°4124, M17SW IRS and VV Ser), or else returned 
values that indicate that although the profile can be a good 
match, the estimated parameters are not physically sensible 
(eg MWC 297). It may be that the success of this model 
in the case of Lk Ha 101 by Hamman and Persson (1989) 
was fortuitous, since it is a very similar type of object. They 
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Figure 7. The profiles of the Fell 1.688^m line in S106 IRS separated by 378 days. The first epoch is the solid line. 



fitted to Ha which will have more easily detectable wings 
for electron scattering, and hence they can be seen to higher 
velocities which is what drives the fit to the electron temper- 
ature. It is worth considering briefly why the model fits at 
all if the parameter estimation is so poor. The reason is sim- 
ply that outlined by Laor (2006) - multiple scattering gives 
rise to a close to exponential profile, and the high velocity 
component of the Castor et al model is essentially very close 
to an exponential. It seems clear that a full Monte-Carlo ap- 
proach to modelling these systems, such as that outlined in 
Kurosawa et al. (2011), is required to proceed much further 
with any analysis. 

We have also shown that the 1.688/xm z 4 Fg-c 4 Fg Fell 
line is seen in massive and intermediate mass young stellar 
objects. We have successfully modelled this line with a sim- 
ple Keplerian disc in all sources bar S106 IRS. The model 
results are largely in agreement with those found from in- 
terferometric near and mid infrared measurements of the 
continuum where such data exist. We therefore ascribe a 
disc origin to the Fell line in most sources, since we know 
from the interferometry and other data that these sources 
do have discs, and our simple model fits are largely success- 
ful. This is in agreement with previous findings from optical 
data that fluorescent Fell emission is a good disc tracer in 
classical Be stars. The counter-example of S106 IRS cau- 
tions against too simple an over-interpretation of this result 
however. 

Overall therefore, these data are strongly suggestive 
that all of the sources observed have evidence that a disc 
is present, as expected for such young sources, and the 
line broadening indicates a very dense circumstellar envi- 
ronment. The likeliest origin for the bulk of the emission 
lines is a combination of a disc and a wide opening angle 
wind or disc-wind. The relative contributions of the latter 
varies widely between different species and even different 
lines of the same species. The similarity of features between 
these sources and lower mass T Tauri stars is notable. 
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